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ABSTRACT 


Measurements of Pu-Be neutron age to indium resonance 
and thermal migration area were made in pure water at 24°C 
pulageliemixtures containing 1.16, 8:96, 27.81, 46.65, 

70.09, and 80.09 grams/liter of boric acid at 40°C. Thermal 
diffusion lengths were calculated from these measurements 
and the dependence on poison concentration observed. It 

was concluded, both from failure of the age measurement 

and the irresolute thermal flux behavior at the heaviest 
Merson concentration, that the "maximum absorption” limit 
was exceeded in this experiment. The absorption cross- 
section of the supporting medium corresponding to the 80.09 
g/l concentration was 0.6120 eae The present results, 
together with previous experimental data, provide an esti- 
Hace sOumciti sion length at the absorption limit of 0.5 cm, 


Pee nvhy less tnan that corresponding to the theoretical 


forneola dimit (0.6 ecm). 
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It has been shown analytically by Corngold [1] that if 
the absorption concentration of a supporting medium is 
increased beyond a certain value (assuming 1/v absorption 
in an infinite medium), no thermal neutron diffusion length 
exists. The practical implication of this is that under 
such extreme circumstances the spatial distribution of 
thermal neutrons cannot be resolved. Table 1 lists theo- 
retical values of the minimum diffusion length which 
corresponds to the "maximum absorption" limit for the more 


common moderators [2]. 
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The dependence cf thermal neutron diffusion length upon 
the absorption concentration in water has been the subject 
of several experiments [3, 4, 5]. In all of these experi- 
ments the limiting value of diffusion length is assumed, 


but not proved. 





The most extensive measurements of diffusion length in 
water poisoned with a 1/v absorber have been those by 
Goddard anda Johnson [3]. Their work included measurements 
in water containing dissolved boron which was ideally 
Biit¢ea for the purpose. Boron is a strong ll Va sonoer . 
i.@,its cross-section varies inversely with neutron velo- 
city, or as the square root of neutron kinetic energy. 
Therefore, the homogenious mixing of small amounts with 
pure water has negligible effect on the average scattering 
cross-section but drastically increases the absorption 
cross-section. 

Goddard and Johnson experimented with a wide range of 
boron concentrations at temperatures of 22.3° and 6526; 
Their most heavily absorbing solution contained 69.3 g/1 
Sue boric acid (H 3B0 2) and gave a diffusion length measure- 
ment of 0.6581 cm (within the Corngold limit). Similar 
measurements using dissolved cadmium at concentrations of 
55.5 and 82.2 g/l resulted in measurements of 0.5983 and 
0,4925 cm respectively, both of which lay outside the 
Corngold limit. They did not report any @hini tae Lille aes sabia 
obtaining these latter diffusion lengths. 

This research project aimed at experimentally determin- 
ing the "maximum absorption" limit by increasing the boron 
concentration in pure water beyond that attained by Goddard 
and Johnson until the diffusion length measurement failed. 

This experiment was conducted using a plutonium-bery lium 
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source which emitted 5.4 x 10° fast neutrons per second. 
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ivi lebe explained in the following sections why this 
Paaencuvmomec@urce OliCtated an unconventional and indirect 
method of determining diffusion length. In the process, 
mitOmmabtion Concerning the dependence of neutron migration 
arpea and age on the poison concentration is obtained for 


the first time. 
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A. TWO-GROUP DIFFUSION THEORY 

(he thermal neutron flux, in most cases, arises from a 
feenrece Of Neutrons having energies many orders of magnitude 
larger than the energy of neutrons in thermal equilibrium 
mernethe medium (about .025 eV at 25°C) [6]. The Pu-Be 
(a,n) source used in this work yields neutrons with an 
average kinetic energy of about 4.2 MeV and the neutron 
spectrum ranges from roughly zero to about 10.6 MeV. The 


calculated oo is shown in Figure 1. [7] 
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- Figure 1. Calculated Neutron Spectrum 
for a Pu-Be Source. 


nie ciiddytacal model which describes the neutron popula- 
tion as two separate groups, one fast and one thermal, gives 


solutions which approximate the behavior of neutrons in real 
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Seenmoltes, ana tor many practical calculations, is suffi- 
ciently accurate [8]. In the two-group treatment, neutrons 
Of thermal energy constitute one group while all those of 
mueiiem energy belong to the fast group. Neutrons are 
eemeraercd removed from the fast group by scattering 
Cormlpetorns which reduce their energy to the thermal region. 
since absorption cross-sections of moderators for high 
energy neutrons are very small, any loss of the fast 
neutrons by absorption is neglected. 

The disappearance of neutrons from the fast group is 
described in terms of a "slowing down" cross section, De 
16]. By denoting the fast neutron flux by on (neutrons/ 


em°-sec), Y. is defined such that 2 (r) gives the number 


F FoR 
of neutrons per em?/sec. which slow down out of the fast 
Geom at the point r cm from the source. 

TiemeimersOveeOp1e DPpOlmy Source at r = 0 emits S fast 


neutrons/sec in an infinite moderator, the fast flux at 


r # 0 can be described by [7] 
Dev bn(r) - Ipb_(r) = 0 (1) 


where v°6(r) is the Laplacian of the flux and Dy is the 
PemicronmeoctiteLent for Che fast neutrons. The source 
ues LoeLaxen into account by means of a source con- 
Timomuiinehe wheel De stated Later in this section. 


The quantity D, can be derived by considering Fick's 


F 
hoaveouwomrtumsion which states that the neutron current 


fnew vemiseaicOport1onal to the negative gradient of the 


el 





fue) Lhe proportionality constant has dimensions in 
memeuneana 1s Called the diffusion coefficient, denoted 
Moeeoeagetne more general Creatment of neutron diffusion. 
Accordingly, the fast neutron current density can be 
represented by [6]. | 

Fast Neut “0 

pee. Density mae DpVop(r) * Fen D(E)V¢(r,E) dé 
where the integral is evaluated over the energy region of 
Miemiascts Troup, 1.e€., from thermal, Ban to source energy, 


E ine fast flux is defined by the integral 


0° 
By 
_ (x) a d(r,E) dE 
E 
ch 
and therefore 
BG 
Jf DGE OCT, ) ds 
Beh 
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J V Oct) di 

Eth 

In two-group theory it is assumed that the variation in the 
House SHeCEruUM 1rom DPOInl tO Point in the medium can 
be neglected, that is the flux is completely separable in 


Space and energy. Thus, 


a2 
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5 ° D(h)o(E)4E 
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Dy, — ce) 


0 
ii o(E) dE 


Ph 


where $(E) is the energy dependent flux. 

Returning to Equation (1), and using the boundary con- 
Sac ton that dn (r) must remain finite as r>~, the general 
solution can be obtained: 


-V/VT 7, 
= & 
on(r) = C r ; (3) 


where the parameter Tin is known as the neutron age to 


thermal energies, defined here as [9] 
fe 3 == (4) 


The constant C is found from the source condition which is 
obtained by drawing a small sphere around the source and 
computing the number of neutrons that pass through its 
surface per second. Considering the magnitude of the cur- 
rent density in fast neutrons per om“/sec, Jn(r), and a 
sphere of radius r, this number is bnr® Jy(r). Then, in 
Peemieaiivwe act => O, Che Source condition is 


nee a Jee) = (5) 
r>0 
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Thus 


do,(r) D,Ce . 
J = - amercmreermrenneerecmommms SS cermemmnnceerpe cee ore 
ian acc ae Fe , 1 
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C= nD, 


ana 


Finally, the fast flux is given by 


~V/ VT 
se 


dp(r) - Wirt, Dp (6) 


Considering now the thermal neutron group, the time- 
independent diffusion equation which describes the thermal 


mis Can be written 
Bye gn(r) - F, op(r) + ag(r) = 0 (7) 


where $,(r) is the thermal .flux, D is the average value of 
the diffusion coefficient in the thermal region, x. iis 

the average macroscopic thermal absorption cross section 

ana An (rv) is the thermal Ss loneins Gown density. The quantity 
Om (r) homegudivatene to the preduction rate of thermal 
neutrons per unit volume and, for the case of a point source 
Oi fast neutrons in an infinite moderator, is equal to 

Enon (r) the number of neutrons slowing down out of the fast 


group. 


Equation(7) can be rewritten as 





Ori 2) meen) 
ge + 2 
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1s 
where dn(r) iomea ven by Equation (2) and Lip ee fies 
a 


wewtenermal neutron diffusion length. By requiring the 








Memtron ilux to remain finite, the homogeneous solution 


becomes 


-r/Ln 


Ky & 


bn (2) y= (9) 


Pees ube source term is proportional to the fast flux, a 


Percicular solution is assumed of the form 


; -r/¥T, 
5 
on(r)p = 2, (10) 


Eyeeoubsetituting this expression into Equation (8), the 


constant Ky is found to be 


-s ae 
Ke, = >: (a 
Am D (Lin = Tn) 


The general solution to Equation (8) can now be written 


-r/L Tey oe 
K, e ae ar a 

bn (Pr) Dine e 30 5 Gi?) 
i Amr D (Lin — Tp) 


where the constant Ky can be determined dy a source con- 
dition as in the fast flux case. However, since the source 
only produces fast neutrons, the required source condition 
for the thermal flux is 


ene eo Jnr) = ies) 
ie 0 
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where Jm(r) is the magnitude of the thermal neutron current 

femsoity. iInus, 

_ dbp (r) 

Jn r) = -p ———— 
Ga. 


and the constant Ky becomes 


O 
7 S Lyn 
“al er no a) 
An D (Lin - Tr) 


Peepoliy, the thermal flux is given by 


7 Bay Ten 
on(r) iia aaa CS -€ ) Gis) 


Equation (15) describes the thermal neutron flux in 
neutrons per em“ /sec in an infinite moderator using two- 
group Peo: 

It follows from the Pe ooeente St eEGueabiom (5) tha 
1t 1S also possible to describe the slowing down of the 
fast neutrons and their subsequent diffusion as thermal 
MeMurons as two distinct processes. The neutron ae. ais 
is related to the distance traveled from the point where 
a fast neutron first appears to the point where it slows 
Sowimuemenermal Cnereies. Similarly, the square of the 
thermal neutron diffusion length, Te COEres ponds to. ene 
distance traveled from the point where a neutron becomes 
thermal to the point where it is finally absorbed. A 


quantity M 7 known as the thermal migration area, takes 


at 
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into account the total neutron travel during moderation 


peemcrtiiusion. Accordingly , 
Mp = Lp + toy (16) 


where Mn (migration length), Tn (slowing down length) and 


La all have dimensions in length. 


B, ONE-GROUP DIFFUSION APPROXIMATION 

Although the two-group theory gives more accurate flux 
solutions, the one-group approximation readily accommodates 
the "first scattering source" theory which is useful for 
our analysis. In one-group diffusion theory all neutrons 
are lumped together into one-group, namely thermal. This 
approximation, while limiting the description of the neutron 
population to spatial dependence, allows attention to be 
focused on the diffusion process. 

In the two-group model of the previous section, neutrons 
entered the thermal group as a result of slowing down out 
of the fast group, and the thermal slowing down density, 
an(r) ; appeared as the source term in the thermal diffusion 
equation. In the one-group treatment, however, when the 
actual neutron source yields fast neutrons, the source term 
used for the diffusion equation must be altered to account 
for the behavior of the fast neutrons prior to reaching 
thermal energies. In the particular case of a point source 
emitting S fast neutrons/sec in an infinite media, the 


slowing down density at r is given by the relation [10] 


dey 





2 
eece JN 
1 (17) 
Ont) - 


where T is the neutron age. A small T means that the neu- 
trons have suffered little slowing down and so they have 
Mermecitiitused far from the source. Larger values of T 
Pemrespond to neutrons which have been slowing down longer 
ena therefore have had an opportunity to diffuse further 
from Bnel sGurce:. When the neutron age to thermal energies, 
Tm > is substituted into Equation (17), the resulting expres- 
sion gives the dispersion of thermal neutrons and can then 
wemised as the source term for the diffusion equation. 
Equation (17) is based on the Fermi age model which 
assumes that the slowing down process is a continuous one 
[9]. In hydrogeneous media, however, the slowing down 
process is not continuous due to the fact that a neutron 
MevelOse all of its energy following a collision with 
Waeecenon sim adaition, the hydrogen scattering cross 
section decreases rapidly with energy above about 10 keV. 
ecm stowecac fast neutrons found far from the source 
are those resulting from high energy neutrons which have 
not suffered a collision. When these neutrons finally do 
collide they will be moderated without moving very far from 
PMemoolmy orecolliston. Accordingly, a fair representation 
of the source of thermal neutrons which can be used as the 
source term in the diffusion equation is the "first 


scattering source" [10] which for a point source is given by 
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a 





q(r) = —, ye 
\ S 
a 


fers equation gives the density of first scattering colli- 
Paemanay distance r from an isotropic point source emitting 
S fast neutrons/sec in an infinite medium. The parameter 
be Pune Macroscopic scattering cross section of the 

mest SOUurce neutrons. The reciprocal of b is known as 

the relaxation length, b, defined as the distance the 
neutron density is reduced by a factor l/e of its initial 
value due to scattering. Replacing Le by 1/b, the source 


term can be written 


K, Pa, © 
q(r) = eo (18) 
r 
where K. = S » aeconstant 
ieee. s. ms 


At large r, therefore, the principal variation in the 
slowing down density arises from a simple exponential 
decrease combined with the inverse square characteristic 
of a point source. Moreover, q(r) is not a function of 


age, since every neutron is presumed to slow down 


Tis 
Pivectemrs first collision. 

The thermal diffusion equation for the one-group approx- 
imation can now be written as 


-r/bd 
o,,(r) K.e 7 
2 L 1 
Voo,z(r) - + ——— —— = 0 (19) 
Al, L 2 O 


T 





Dr 


ny 





To satisfy the boundary condition that On(r) go to zero as 
r+ ©, the homogeneous solution becomes 


on (r) 4; = aa (20) 


A particular solution is then ASosumecdso! che Ler 


ee” ae" 
dnl) p a ey + an. aan ° 2A.) 


J 


By substituting this expression into Equation ClO Re ene 


constants C, and Co are found to be 


and 


a (22) 


The general solution to Equation (19) can now be written 


-r/Lp, K “ai 2~r/b 


Ke il ae 
DPA Gly ee (23) 





fartic  adurtuswen length, Ln > foelercer than the relaxa— 
tion length, b, the first term of Equation (23) dominates 
at sufficiently large r, the thermal Tluxeeecays expo- 
nentially, and Ln can be evaluated in a simple routine 
manner. This is the conventional method of determining Ln 

The relaxation length, however, depends on the initial 


energies of the neutrons generated by the source. Fast 
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Peuwerons from a Pu-Be source have relaxation lengths of 
about 10 cm. in water [11] which exceed 2.78 cm., the cor- 
responding diffusion length of thermal neutrons [6]. Con- 
squently, the thermal flux decays in a more complicated way 
and the simple theory for the case of a completely thermal- 
ized neutron field does not apply. The JeveriiEnavlom oO i Lr 
iimwaver can, however, still be accomplished by measuring 
Pe@aethne Migration area and the neutron age. From Equation 
(16) the diffusion length is then obtained from the relation 


by e r 


T t) ou 


CC. DETERMINATION OF NEUTRON AGE 

The neutron age, tT, between two energies is proportional 
to the mean square distance travelled by a neutron as it is 
being slowed down from the higher to the lower energy. 
Denoting this distance by r as the second spatial moment of 


S 
the slowing down density can be written as [9] 


_> f r°q(r,1(Ep,E)) av 

Po (E, >E) = Far Ey BE) av = 6T(E,,E) (25) 
where Eo is the energy of the neutron at its point of birth 
as a fast neutron and E is a lower energy attained as a 
result of the slowing down process; dv is the volume element, 
and q is given by Equation (17). Neutron age is now defined 


by the relation [9] 


Se 
2. 


T(E, 5E) = é 


(Ep sE) (26) 


Equation (25) applies to neutrons of any energy and, 


hence, of any age. In actual measurements, however, the 





Pirie ace is 11xed by the resonance energy of the detect- 
ing foil. Consequently, in determining the age to thermal 
energies, Tm > the spatial distribution of the slowing down 
density at some energy above thermal is obtained and then a 
Correction is applied to account for the additional age 
megured CO reach the thermal region. 

The slowing down density of source neutrons at 1.45 eV 
Can be determined by activating cadmium—-covered indium foils 
peecea at various distances from the source. The principal 
indium isotope, ie hoiomomoticonees 1SsOlaved absorption 
resonance at 1.45 eV while the absorption cross-section of 
cadmium is large below 1.45 eV, but small at higher energies. 
The cadmium cover, therefore, prevents thermal neutrons 
from reaching the indium foil; whereas the large indium 
resonance at 1.45 eV is primarily responsible for any 
Meutron induced activation, 

The saturation activity of a cadmium covered indium 


ime 


foil, A iomedUal yOwuhe rave Of production of IN 


CD 
nuclei and can be written as [9] 


Aap = Vf o(r,E) Fay dE (27) 


where V is the volume of the foil, de menunhe macroscopic 


IN 
oe ands the integral is 


absorption cross-section of IN 
evaluated across the 1,45 eV resonance. The number of 
heutrcencm per om?/sec WilLCheestlewraQowh past a given energy , 
tions LOMm approximation, is given by 

q(r,E) = $(r,E)& 2 (E)E (28) 


where $(r,E) is the energy dependent neutron flux, i CE) 


is the macroscopic scattering cross-section, and € is the 


oS) 





average logarithmic energy decrement. From Equation (28) 


meer Oollows that 


q(r ,E) | dE 


Aap = V a F (29) 


Pesuming that q(r,E), €, and z (E) are constant over the 
eel Cnergy range of the resonance peak, this can be 
written 


A _ We Gea. f r 


aE 
OE) == (30) 
CD Er. ‘ E 


i, 
That is, the slowing down density at 1.45 eV is proportional 
memune experimentally measured activity of the foil. By 
substituting this result into Equation (25), and noting 

that dv = War-ar, the expression for the age to indium 


meoemance 2n Spherical coordinates is 


oa zB Hi 1G Aan dr 
T = S = te Oe (31) 
IN 6 6 i 5 
Sf or Aap or 
0 


Thus, by knowing only the relative values of q(r, 1.45), 
the age to 1.45 eV can be obtained by numerically evaluating 
the integrals in Equation (31). In this analysis the slowing 
down density is assumed to be approximated by Equation (18). 

Another method for measuring neutron age [6], which 
provides only a rough estimate in hydrogeneous media, is 


Pa wmoamiien spatial distribution of neutrons of a given 
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em@erey described by Equation (17). Taking the logarithms 
of this equation and considering indium resonance energy 


leads to 
- o 
Qn q(r) = const - r /Atry . (32) 


Since the resonance activation is ee ional TOncim) a 
gmot Of &n Aon measured at various points from the source 
can be made. For moderators which satisfy the continuous 
Slowing down theory of the Fermi Age model (e.g. graphite) 
meacmexperimental points will lie on a straight line. For 
water, however, the points will lie above a straight line 
excep. in a very narrow band about 20 em. from the source 
Piel. this is due to the unusual nature of the hydrogen 
secattering cross-section and the fact that individual 
neutron energy ligcecemeweusOmlacee as previously discussed. 
iMmevars narrow band, however, an approximate straight line 
can be drawn whose negative slope will provide an estimate 


OL Try 


D, DETERMINATION OF THERMAL MIGRATION AREA 

The thermal migration area, defined by Equation (16), 
correlates with the total mean square distance a neutron 
Maewetsmironedts DOlNt Of birth to its point of capture. 


This is given by [11] 


rr = (33) 
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where Ag is the total saturated activity measured as a 
function of distance r from a point source. Equation (33) 
Gan be obtained in a manner analogous to the derivation 
Smeeguation (31) for PS Here, however, it is the total 
memarave neutron flux which must be known to evaluate the 
integrals. At energies that activate a bare indium foil, 
mere boval flux will be proportional to the saturation 
femvicy. Migration area can then be determined from the 


definition : 


a 
Mio = (34) 


ae 





Peete ON OF THE EAPERIMENT 


Measurements were made in pure HO at 24°C and in 


O 
imerures comvainimg 1.16, 8.96, 27.81, 46.65, 70.09 and 


SOe09 grams of H BO. Beo@mibittemeomweH 0 at 410°C. The wast 


3 2 
TWO Mixtures contain the highest boron concentrations 
movained Co date. 

Mie “erinieneosweren comaucted in a cylindrical plexi- 
ess bank, 1.45 ft. in diameter and 4.27 ft. high, which 
held about 50 gallons of water. This tank was situated 
puopearsmMalléer, steel tank 2.95 f{t. in diameter and 1.25 
ft. high, which held about 65 gallons of water. The base 
tank served to contain the Pu-Be source and also provide a 
protective water shield. When positioned for measurements 
in the base tank, the source was 0.625 inches below the 
Pos wou uieomolex ie lass Gank at centerline. A polyethylene 
weeneupeCnwleconin. thick made the source holder watertight. 

A single steel rod immersed in the tank at the side 
heated the poisoned mixtures through a rheostat control. 
Uersmicaminemroa Olas 4 Continuous running electric stirrer, 
maintained the medium at a uniform temperature of 104°F 
Es ae 

The relative flux was measured with bare indium foils, 
IMOmineiineciameter ang .0O10 in. thick. Foils used in any 
given measurement had a weight spread of less than 2b. 


Indium resonance activity was measured with the same foils 
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milcced Inecacmium covers 0.025 in. thick. A plexiglass 
foil holder arrangement provided accurate spacing of the 
fools in ewehesmeatum On a line perpendicular to the source. 
The distance between adjacent foils was 0.50 inches for 
mimeimeasurements. Figure 2 shows the experimental appar- 
atus and set-up. 

Nine foil positions wereused for each measurement. 

The closest measuring position to the source was 8.26 ecm. 
while that furthest away was 18.41 em. Counts of the cad- 
mium covered foils became insignificant beyond this last 
Pes tLon. 

The cadmium-covered and bare indium foils were irradi- 
ated separately during each experiment for a minimum of 8 
hours. After irradiation the indium was allowed to decay 
ieee 3 minutes to avoid interference with the desired 
ie lccwaie Grom the 13)second half lire of in1?° 
[11]. The activity of the foils was then counted from each 
side Tor Memimnuves: the relatively long counting time 
having been chosen to improve statistics. The counting 
system consisted of an end-window Geiger Mueller tube con- 
tained in a lead shielded housing, and an eiectronic 
Scaling Unmet sow apparatus 1S Summarized in Table 2 
ands piLeCturcdmaum Figure 35. 

Counter background was measured repeatedly during each 
run as an erratic line voltage perturbation appeared at 
infrequent intervals. Only rarely was the background 


Spable within statistics for an entire run. The resulting 


of 
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xperimental Apparatus 


Figure 2 
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TABLE 2 


DE TEC een Eh OULPMENT 


Counter 
Activity Detected Beta and Gamma 
Type Organically-quenched GM 
shape Cy aaacies Cra) 
Counting Gas Helium 
Manufacturer Atomic Accessories, Inc. 
Polit Ter 
Type Non-overloading 
pensitivity 10 mV 
HeanuiacLurer Hew Lemusbacktard 
Model ‘ 5554A 
Secaler-Timer 
ManinacLurer Hewlett Packard 
Model S201 


necessity for interpolating between background determina- 
CLONS represents one important possibility for nonrandom 
Eauors im the results. 

Prior to each poison experiment samples of the solution 
Veuomencmile- Livesanalyzed for boric acid content. The method 
Specdalyalomvas DOrdce acid titration with glycerine which 


is outlined in Appendix A. 


30 





ive ena otS OF DATA 


A. DETERMINATION OF SATURATED ACTIVITY 

Three corrections were applied to the foil data. The 
Meesb Was an internormalization based on foil weight, the 
foerewo COrrecued [Or absorption of indium resonance neutrons 
by the cadmium cover; and the third corrected for counter 
resolving time. 

The weights of the indium foils used in all measurements 
varied from 0.9282 to 0.9621 grams. All foil weights were 
therefore normalized to a weight of 1.0 gram and the result- 
ieecOrrecrions applied to the foil counts. The highest 
Value of the weight correction factors obtained in this 
manner was 1.082. 

The indium resonance counting-rate correction was 
fear wine a to be |i. Ins value was obvained from a 
wie eOmscorrecti1on factor vs. indium foil thickness con- 
Paticreadsnor a cadmium thickness of 0.025 in. as outlined 
Puetitule (14). 

Mm number of determinations of counter resolving times 
were made by following the decay of 54 minute tnit6 over 
three hours. An average value of 400 usec. was obtained by 


solving 
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where A is the extrapolated activity, N is the mean counting 
meen and t 1s the resolving time. A correction factor of 
1/(1-Nt) was then applied to all counting runs. 

The effect of foil activation by source neutrons above 
resonance was considered negligible as the closest measure- 
ment point was greater than 8 cm. or the source [12]. 
Mmieoreco insure that the close spacing between foils did 
not contribute error due to shadowing, activations were 
made in pure water using various spacing intervals. 

VaMeincomrecolie erie 1Oll count as described above, 
the saturated activity in counts per minute was calculated 


heom the equation 


os AC 
~ - =_ = — — 
S At] =A(E-b,) ACE 3-t) 


where C is the total count less background, A is the disin- 
tegration constant of rnil6 ty is the irradiation time, 

Co ty is the waiting time and eal is the waiting time 
plus counting time.[13]. An average quantity was then 


obtained from the values of A, computed for each foil side. 


5 
The computer program used in this calculation is included 
in Appendix F. Appendix I provides tables of the computed 


Gata for all runs. 


B. DETERMINATION OF RELAXATION LENGTH 
It was pointed out in Section II-B that in hydrogeneous 
media, the slowing down density at large distances from the 


source can be approximated by the expression 
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reo ae 
a(r) = =—5— (18) 


ie 
where b is the relaxation length of the fast source neutrons. 
Since in the framework of age-diffusion theory, the flux 
is proportional to the slowing down density at a fixed 
energy, it follows that the cadmium-covered activity at 
large distances from the source would have the same space 
dependence described by Equation (18), namely 


= 
Agp(r) = +4— (35) 
Ir 


To evaluate b from the corrected data, the product 
Aap? was plotted against r on semi-log paper. This plot 
showed the spherical Space distribution beginning to 
decrease exponentially at about r= 12.5 cm. The distribu- 
EEo@eLor 2 > 2.5 cm was then fitted, by least squares, 
with functions of the form of Equation (35); specifically 


Qn Aapr® = const - r/b (36) 


Hiemousi@tste tome tor r > le.5 em consisted of only five data 
Pebaecme whiten, tor most runs, exhibited large statistical 
Vaularnonmawe primarily to the low foil activity in this 
range. Accordingly, the solution resulting from fitting 
(36) to this data provided only a rough value of the param- 
eter b. A more accurate determination was made by comparing 


Aap distributions calculated from Equation (35) (using 
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neighboring values of b) with the measured distribution on 


a large scale semi-log plot of A VseeeGes Unlike the 


CD 


‘oa spherical space distribution all data on the Aon 


Mmeeoun decreased in value with increasing r; this permitted 


Vict 


Meemworetme more reliable data close to the source to indi- 
@aee which, if any, of the data at r > 12.5 cm were incon- 
Sistent with the decreasing trend of the overall distribu- 
bon. Atter discarding any poor data indicated by this 
procedure, the calculated Aan distributions were normalized 
and fitted to the measured distribution for r > 12.5 ecm 
Mitel one Closest possible fit was attained. The value of 
Peo iim Calculating the distribution found to give the 
elosest fit became, then, the negative reciprocal of the 

re Visnsorcurve., 


or 


The above procedure was also used to determine relax- 


slope of the exponential portion of theA 


ation length for the bare foil activity (A,) and, thus, the 
Sapesot the Agr® metic mm Ostet. lm measurane® the 
age of Na-Be neutrons [15], assumed no loss in precision 


from using the same value of b for both the Ag and Ag 


D 
Prispripyutons, The present procedure, however, clearly 
showed a small difference in the two, the b for the 
resonance distribution being less in every case, as it 


should, The work of Rush with Ra-Be neutrons [16] showed 


CD cS oenie ea 1G) 


Ag VeenOor sine pure Water run which includes the data 


this same tendency. Figure 4 is a plot of A 


computed from Equation (35), normalized to the experimental 


data. Figure 4 also shows the marked deviation between 


34 





TUM rT ery 


Te 
eth thay eT tr 
Te eee 
| net pepe ae aE ee 
eee wf 
i Ce et 
i ee eee i 


TTT [An el 





11] at 
ze | 
THT th) 
‘ | ee 
| ) ) 
de 





' - : t . , { eee = 
i bi = i 
—4{ + - - 
+ 





Sf ia t 
ra oi : : 7 

. * e 
acd DUT LO! 


ean tact oe 
G 
4 ‘ t + . 5S 
| : a i . e+ 
; 
si 3 
= + 
+ 





fee ree 


SS ee , 








ae a ae 


ce erect ee : 

na oC TEE aa HHI te ee Ha 

Poe aan br HE Pee 

eats tt CoA 7 ) mT a 
pe CECH ic Ht 

rH Sa ari 

ete oie Heer A 

Soto EOE EE ee Me Lh 

; POE COTTE CETTE CECT 

ieeemieeee A =e Heer 





AT tite SO CC CTT qt 
Ree At oo 2 SEY meas Pete tt Hirt rT 
TET iu Sap SCOTT po | SST a 
RAT eo a eT ACT oT 
Pon yt a a TTL PARP ee 
aAREt Chy/-4 OMAN TSELTOR DOCS EERIE Bem OR On HTT LSE CE co re 
SAA ro cn Eran a Tee eT 
Se ttt ret CO 
Mier g et ae me UL TTT PC SET 
aye Laika EM ALE STEPH Py 1 leas a Sa co 
os (do) IAG Ea ON aI ae 


- cl 


os 
Sp 


Figure 4 





the computed and experimental data for bermars mear the 
source. This is to be expected since Equation (35) applies 
only at larger distances from the source. Similar plots 
for the poisoned water runs are eee in Appendix B, 
The values of b determined for all experimental runs are 


Paevem in Table 3. 


TABLE 3 


RELAXATION LENGTH 


02030 cm) 

Berne acid Indium Resonance b(+£0.02 cm) 
(B72) Activation Total Activation 
--- LO2 50 10.50 

eo + . 17 POS 5 10.45 
8.96 +°. 36 lf On.25 Orns 6 

Poa ae | 9 Oren 10.36 

HG265 f2h2 TOUS nO 5.35) 

7fO.O09rs 1.3 HOS Or on) 

80.09 +1.9 EO eles OSS 


_—O > rh 


C. DETERMINATION OF AGE AND MIGRATION AREA 

A smoothed curve was drawn through data points less than 
r= 12.5 cm. on the Neo and Agr® Weigabiss i jolene. | lene 
rb 12.5 cm. the data were extrapolated to infinity by means 
of a straight line having a slope equal to the negative 


reciprocal of the parameter b. The areas under the smoothed 
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mm veo were found by integrating the straight-line portions 
mm@alyoically and by using Simpson's Rule for the remainder. 


Points from these curves were used to obtain the curves 


4 


Aopt and Aer’ Versuowiemwinese Pacver curves were found 
to approach straight lines beyond r = 45 ecm. and areas were 
Careultated accordingly. After the area under each of the 


four curves had been determined, the second moments r # and 


= 
PS were obtained from Equations (25) and (33). The neutron 
age to indium resonance and total migration area then 


followed from Equations (31) and (34). The computer program 


ieedeith Ghese computations is included in Appendix G. 


2 2 y y 
cpr > Ag? 3 Aap? weand Ag? 


Memecus tf fOr pure water. Piots of Acpr® and Agr* versus 


r for the poisoned experiments are given in Appendix C. 


Figures 5 and 6 are plots of A 


Table 4 gives the values of Try and Man determined for 
all experiments. The errors cited have been estimated from 


a and Aar® versus r plots by using the individual, 


unaveraged foil points to construct maximum and minimum 


2 
IN and Mn elle 


A rough check on the experimentally determined values 


values for T 


of age was made by fotlowing the procedure discussed in 
Section II-D with Equation (32). Values of activity as a 
funetion of r had been previously computed from Equation 
(35) in determining the closest value of relaxation length 
to fit the experimental points. This data was plotted on 
semilog paper against r*. A narrow band between r = 18 


and r = 22 cm. presented the closest straight line 
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TABLE 4 


MIGRATION AREA AND NEUTRON AGE 





ES 





Boric Acid 2 2 e : (om® 
g/L aoe "inten ) Fee a, =, 
--- Ciao? 1, cle 52.19 40.79 Se 5 
Mero 17 86 459.41 +0.68 52.064+2.34 pam 
ooo + , 36 boo +1. ly emus. 16 pile oe 
Beer = .79 54.2142.89 51.96 + 3.06 pi 55 
Mes .2 53.47 43.13 51.29+3.42 © 51.47 
moog +1,3 Been ee les Sileaciont. os il lee oe 
80.09 +1.9 53.13+43.13 51.20+4,21 Duleesis 


eee 


approximation. The slope of this line, equal ton i 


iE 
‘from Equation (32), was then determined by the Wie ct 
east squares. <A plot of Aon vs. re for pure water with 
b = 10.3 is provided in Appendix D. Table 4 includes values 
of Try obtained in this manner for all experimental runs. 


The least square machine program used in this determination 


Steaese wis ancluded in Appendix H. 


D. DETERMINATION OF DIFFUSION LENGTH 

Prior to determining Ln» from Equation (24), it was 
necessary to apply @ correction to Try to account for the 
additional age required to reach the thermal energy of 


0 .025i-eV. 


40 


ss 6@ P= 





Pawiiiaves Ol Lhis Correction for hydrogeneous media are 


2 


found in the literature ranging form 0.3 em” [17] to 1.8 em 
fii}, the latter value being that used by Valenti for 
Teaeurements with a Pu-Be neutron source. In this investi- 
gation it was found that adding a correction of 1.8 om* to 
the value of Try determined for pure Cee yielded a value 
of 2.786 for Ly» This is in excellent agreement with 
published results obtained by the more conventional method 
Giemeasuring the exponential attenuation of thermal neutrons. 
Mmeeas also noted that 2.786 cm was in precise agreement 


With the value of Ln obtained from the empirical equation 


[14] 


Lm = 2,64 + 0.0061 T em (37) 
H50 


inerc a temperature, T, of 24°C. Accordingly, a correction 
onl, om” imvuic ace. | TOneimnilumeresonanice LO CThermal ener— 
Giese Was applied to the experimentally determined value of 
ty, for all runs and L, was then computed from Equation eel 
Results are presented in Table 5. Errors have been estimated 
from a combination of those previously cited in Table 4 for 

2 


c 
IN and Mp . 


E. DETERMINATION OF THE THERMAL FLUX DISTRIBUTION ACCORDING 
TO TWO-GROUP THEORY 


Peiigeemeneck on the results obtained for pure H0 and 
the poison mixtures was made by computing the thermal flux 
distribution from Equation (15). The experimentally deter- 


mined values of Lin and Tn were used in the equation and the 


4} 





TABLE 5 


MEASURED DIFFUSION LENGTH AT 24°C IN 
PURE WATER AND 40°C IN POISONED WATER 


—————————— Eo 








Bowler cid 


e/2 Lin (cm) 
0 ee ae 

ae | 2.7862 + .14 

leetcmeceel Ae Oe 

Cecio 30 Ase eee 
Piero: -. 19 Ono OO4 29 

46.65 + 1.2 G.6203 2 .20 

Or Olo mesa des 5 0.5594 + .18 

SUn 0 Omse I. 9 Oe oe 5, 30 


RE 


total macroscopic absorption cross-section of the poisoned 


medium was computed from 


Da = 2. ee 
hydrogen HOROn 
where ye was taken to be 
hy aroren 
» = 2 + N x 30 
ie Be 
“ny drogen eEao Se “nydrogen 
and 
y = N ¢ © (0) 
2Boron HBO. 2B oron 


The parameter o, is the microscopic absorption cross-section 


a 


and N is the number of molecules of HBO. per em> 


HBO. 


Getermined from 
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A 


eo. 0 


. q Burs 
Fes) ; 
ey Molecular weight 


where op is the density in erams/cm> and Ag is known as 


Avogadro's number, which is 6.024 x 10° molecules per gram 
mole. 

The computed values of flux were then normalized to the 
thermal saturated activity which was obtained by subtracting 
iiemeovered foil from the bare foil data. The results were 
meee ueca Cogether on semilog paper and compared. Figure 7 
mies this comparison plot for pure oe Sulake em je lene 
for other runs are contained in Appendix E. Table 6 lists 
the values of z determined for each run together with 
the maximum deviations found between the experimental data 


Sae@meiix distributions obtained from Equation (15). 


TABLE 6 


COMPARISON OF EXPERIMENTAL RESULTS WITH SOLUTIONS 
OF THE TWO-GROUP THERMAL FLUX EQUATION 


a. 5 (om™*) MAX deviation 
--- OrmiOze 2 la ey 
eowe. v7 C030 7 age 2 
8.96 + . 36 0.0882 eS 
Oi ol, -£ack7 9 032270 ae 
NG) 6 Ge ae al Os oo 5a 
Ono +. i. 3 0.5384 5 
SO Ome a9 0.6120 --- 


Nene ee ee meee ee renee eee ere een nnn een enn nn ee ere eee cr eee ener en ne a NN Se 
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ieecOMeanmicoON OF RESULTS WITH THEORY AND 
PREVIOUS EXPERIMENTAL MEASUREMENTS 
A. AGE AND THERMAL MIGRATION AREA 

Theoretical investigation of the age of fission neutrons 
in pure H50 has been extensively pursued and, for this case 
pmeene. there exist calculations of the actual slowing down 
distribution which can be compared in detail with measured 
Gistributions. Comparable results have not been obtained 
BeieMor—11SSi0on sources due, primarily, to uncertainties in 
Miemenerety spectrums. There 18 no account in the literature 
of an age calculation for a Pu-Be source. Coveyou and 
ela vane (17), however, have calculated the age to indium 
resonance, Try,» for monoenergetic sources in H50 using a 
Monte Carlo age code which provides a rough comparison for 
mibeeswork, The plotted results of their calculations demon- 
strate a rapid rise in age with increasing source energy; 
and for an initial neutron energy corresponding to the 
average energy of a Pu-Be source (4.2 MeV), these results 
give a value for Try of 52 em* 5 em). 

Although clearly defined, thermal migration area, Mn 
Pelee iend re miomca McWlalion because Of theoretical 
difficulties encountered in describing the scattering of 
thermal vreutrons [18]. As was borne out by this experiment, 
however, thermal diffusion length is small in highly 
absorbing systems so that the predominant contribution to 


Mao comes from the slowing down process. This allows some 


us 





comparison of measurements of Man aicescaleulations Of Tm 
mee reaving tre thermal diffusion as an additive correction. 
The only measurements of Man and Try from a Pu-Be source 
which have been published to date were those by Valenti 
eae oullivan for Try [19], and by Valenti for My [11]. 
mhneir results are compared with those obtained in this work 
metaole 7. A factor which accounts for some of the dif- 
ference in these two measurements is a source effect cor- 
rection used by Valenti and Sullivan [19]. This involved 
meerererectton to the tndium resonance activity to account 
Miee@eniinite Size of the source and was intended to 
Mecenimave a GCistribution closer to that expected from 
a point source. Due to the relatively low counting rates 
encountered in this work, it is estimated that a source 


effect correction would have increased the activity of the 


foil nearest the source by less tetarcr ieee ll) 


TABLE 7 


COMPARISON OF RESULTS WITH PREVIOUS EXPERIMENTS 


Plutonium—Beryllium Neutrons in Pure Water 





Valenti & 
Present Results Slley ata.) 19] Valenti [11] 
2 2 
Mm (cm) Gules forte lee G4 61037 
Try (em™) 52.19 + 0.79 SOOO 
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In Figure 8, measured values of the parameters Man“ and 
Try are PUGtread 2582 'Unet1on Of poison concentration. The 
Siiterences bevween the two curves is seen to reflect the 
1/V behavior of the boron absorption cross-section. At low 
Thermal energies absorption is heavy and affects only Mm : 
femcrerfeies above thermal, absorption falls off rapidly and 


has a very small affect on both Mn” ST 


B. THERMAL DIFFUSION LENGTH 

The close agreement between the theoretical value of Lin 
mr pure H0 at room temperature with that obtained in this 
experiment was pointed out in Section IV-D. 

Previously reported investigations of Lim dependence on 
maucmpeison concentration of H0 were conducted using a ther- 
mal reactor as a neutron source [3, 4] and a Sb-Be (25keV) 
source [5]. As a result there existed only a small com- 
ponent of fast-flux in the region of measurement which 
could easily be eliminated from the experimental data. 

This permitted relatively straightforward measurements of 
the exponential attenuation of thermal neutrons as discussed 
Mmecerrtonel hCG and, consequently, gave final results of 
eee accuracy . 

Pigure 9ecompares the results obtained by Goddard and 
Johnson [3] with values of Lp found in this experiment. The 
two sets of data are found to be in good agreement when the 
pee temperature difference is taken into account. Since 


their results were obtained at 65°C, their values of Lan 
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mioulompe Migher e@ecording to Equation (37). It can be seen 
from the curve representing the present measurements that 
Boor bing solutions with concentrations of 70.09 and 80.09 


ey ieilay outside the Corngold limit (0.6 cm). 


C. TWO-GROUP THERMAL FLUX 

The thermal flux computed from Equation (15) is compared 
with measured values of saturated activity for each experi- 
mental run in Figure 7 and Appendix E. Due to the compli- 
cated manner in which Pu-Be neutrons slow down, these 
Pempuvcauions provided the only method of examining initial 
activity data for agreement with a distribution predicted 
Peeermcory, tf 1s evident from the maximum deviations given 
in Table 6 that the measured and calculated distributions 
come closer. in agreement with increasing poison concentration. 

The purpose of the experimental set up was to approxi- 
mate as closely as possible the "infinite medium" condition 
Wikeie wag a Dasic assumption in the theoretical treatment 
of both age and diffusion length. That there was, in fact, 
no leakage in the axial direction is deemed a valid conclu- 
sion from the absence of detectible bare foil activity 
bevyonasso.cm. trom the source. An indication of the extent 
of lateral leakage was obtained by measuring total activity 
CaieeyNouULrwarG from vie Source in pure HO. The extrapo- 
lated distance beyond the physical dimensions at which the 
flux can be assumed to fall to zero was found to be 1.29 cm. 


Pieter coulLvmis GCloeplayed in Figure 10 from which it can be 
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seen that even in pure water the "infinite medium" assump- 
TiOnN was not seriously violated. The minor differences 
between the measured and calculated flux distributions for 
eoncentrations of 46.65 and 70.09 g/l indicate the better- 
ment of the assumption at the heavier absorber Concenvrarion. 
toe tsemost apparent from a study of the plots in Appendix 
Ere 


ce 








VI. CONCLUSIONS 


ine dependence of diffusion length on poison concentra- 
mien Was found to agree well with that reported by Goddard 
and Johnson [3]. However, because diffusion length as 
eorained 1n this experiment is the square root of a small 
difference between two large numbers, the results given in 
mabe > Ineclude Sizable estimates of error. Errors in the 
measurements from which diffusion length was obtained were 
due to limited experimental conditions, the most serious of 
Pieehe was the weak neutron source, resulting in a low activity 
am@oea small number of data points. This Pie turhie Miereases 
Mien sbavistical error. 

Previous measurements of neutron age and migration area 
from a Pu-Be source in pure water [11, 19] agree with the 
present results within the experimental errors. This report 
is the first of measurements of these parameters as a 
Miemcimoumsoorsom COoneenuraty;on. Accordingly, the results 
given in Table 4 should be of considerable value in extend- 
ing our understanding of the slowing down process. The 
Peto mass eerabedewivh cach measurement Nave been conserv- 
atively estimated and are meant to serve as bounds on the 
Suet rumici = tiaam a probable error. 

It is evident from the experimental data plotted for the 
heaviest concentration (80 g/l) that it was impossible to 


make a reliable age measurement in this case. The value 


ZS 





mavens tOr this run is therefore merely a logical estimate 
based heavily on the previous measurements. It follows from 
Moor iat vbne numerical value listed for diffusion length 

au this concentration is also unreliable. However, since 

it can be assumed that each of these values approximates 

the result one would expect to obtain had a measurement been 
possible, they were used in Equation (15) to calculate a 
mieieadetlux agistribution. The result of this calculation 
imeloucea in Higure 1] together with the experimental data 
peroeUeh Which no curve could be drawn. Figure 11 vividly 
shows the irresolute behavior of the thermal distribution 
memene oo 2/71 boron concentration. It is concluded from 
Veena ttom the failure of the age measurement at this 
concentration, that the "maximum absorption" limit was 
exceeded Tree experiment. The cadmium measurements of 
Goddard and Johnson together with the present results 
iaemeate that the diffusion length corresponding to this 
miei momapeur 0,5 cm, slightly less than the theoretical 
value, The absorption cross-section of the supporting 
medium at the time the final measurements were attempted 


was 0.6120 a 
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APPENDIX A 
METHOD OF CHEMICAL ANALYSIS FOR DETERMINING 
BORIC ACID CONCENTRATION 

A. PREPARATION OF STANDARD 

A legram sample of HBO. (Mol wt 61.84), finely ground, 
was dissolved in 300 ml of pure HO. Meme iisSe solution ana 
beaker were added 10 drops of phenolphthalein indicator and 
25 ml of neutral glycerol. From a buret, standard NAOH was 
added dropwise until a distinct reddish-pink color appeared. 


This marked the true end point for boric acid and established 


the molarity (Z) of the standard in moles/liter. 


PeeeeaANAiyols OF MODERATOR CONCENTRATION 
To 25 ml of moderator solution were added indicator and 
NAOH as described above. The standard alkali was then 


aCe edeopwise Unell the true end point was reached. 


C. EXAMPLE 


imeomioreiwnOn preduced end point in solution 
containing 0.162 moles HBO... 
7,001.35 


Ofs ml of standard alkali produced end point in 
25 ml of sample drawn from moderator 


LS = 12.15 moles/liter standard alkali 


(25 ml)(molarity of sample) = (.3 m1)(Z) 
molarity of sample = .146 moles/liter 


Cole od) (2146) = 9.03 g/l HBO. Coneeneravidon of 
moderator 
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APPENDIX D 


PLOT OF CALCULATED INDIUM RESONANCE ACTIVITY 
VERSUS r2 FOR RELAXATION LENGTH OF 10.3 ecm 
































ae 4 ~ 


+ - pals = 
eee eed eee - 








































































































2000 


1000 


hOO 


Figure 24 


LO 


5 
ee 





APPENDIX E 


PLO sOr thE nite mtx Vou FOR POLSON RUNS 

















fe ISe ean - me eee ee fae ts 
eee 
* 
— r 
+ - I—- ~ -++ - 
+ 


ae 





rma 











A-+— Wa 
nie 


" asure 











ef > 


























| a Pore | pets 
| i aee a 
PE a Eats . cH 
a tf sa tanaseepfeaseae 
SE a reer ER SSPE 
ti i it aT a Oe 4 bee eeestneeee Te eeagar 
Taree yt RO 








Hitt uTT SC ieee 
cara ea eae cet 


Ht mec ECoe 
| Soiree ea ae 
0 n N 





ia 


10 


r(cm) 


Figure 25 


70 








plrabate 
{ \ ! 
my Hitt ona 


1 deg thor ed 4H eM eat ih lies tu Leal 
reterresterttrerrerestees tippy tran crea 
{ ta Tice a {| cals | | 
riptrriory tetesteterf ett | : 4 | | | 
+ | F ' ti ‘i Mi if ' | 


r ' 
le 4 
‘ 


4 re 








| 
































+ sae Ne 
A - } 
4 cabs 
= = + 
ia s * 
— } 4 
= aN 


r(cm) 


























“ica 








FETE 
4 FEE 
aE 





a 


vw 


Ht 
TES EEEEEE PHS 


He 


| 


ae 
aaiicea 


© 


| mt | 
Se etait | saa 
ie 
a 


Be ee 
CEE 


SSE ISS 


ON 


= 
co 





10 


Figure 26 


ffl 





















































TT pty 
! Hitt 
iT B 
n fh HH 
tel 
Sut 
; NEG Hy 
TTT 
eh 
OO Tt 
MME po 
TS ih ee aE 
TT eT Heth | 
i esirt int 
4 ths i ++ 
mic: 
ity a 
re 
: th 
© 
Ot: . 
{ Lt 2 
ROOT bi S 
TTT a i 
fr, 
oO. 
temo 
af 2 
mm jo 
+s 
a nT SIsUEUERESRSSSTOGGEE EGE 3 
mm ey tt Daa ele ee Coty Bis 
iN Se Tr r ASh ed 
eit [peat + | een 
i fot Th itSRESuu PSE fo 
TT ac Po rE cot ae | ce oe 
| ae ut BH rah Etats i ers 
Enh Pe oe Pata 
rye Leth f a SaaS ee a at alee ect cana See 
“ te | ye fe. let eed Pett h ed Lt ete 
Cm MTT TE REE PE fz tan 
2 CCT eTEEEEEer TP 
Tp Ay Ty Mina Poorer PTT ferry y op 
0 6 |6O © K © yt = i 0 « 
0 i 5 
<[ <7 





ie 








coo r Ta ESE ; 























r(em) 


Figure 28 


73 


























gece etree 
| EES Ee 


N 











“Og 








SEE 
ty —: 


‘Aetivi 


a 
nae 


{id sO. ut! 
ration-70..09 





= 
i 


a 
ae a 
C cent 


Hlty 
oe 


+f 


cent 


ae) 


ze 
- “ 

Therma 
St 
onc 


eee 
cea tnee taal ek 
{ He z 


ee et 


Ge 


= le 


nc 
easure 


et 
= es 


a 
2PM 


























ae 
} 


awe Cy 
' . ’ 


ty 
' 


vprert Pett fors 
bere deren 


rae the fy ate Wd 

eae ts i 
{'i'aeepaf || Pad 
ct 

oth ttt \ 


TR TT 


t at 7 
ieee clne if! i 
ed, eh ies ' 
| tle lt: ; | a rl 
r w—erh tp : 4 44 + 
"ue debe to rode i j ! 
I Weeetadhs: aad, barre pl i 
» + '" ' 
a : ! at 
- 
rf : 
' 


tet 





ye tla 
batt ated ied fa 
‘ atte aaonhe 
| apts ts ' 


, | 
11 

sport per epeees 

' ' 
fs bit t , 

Neate ty 
Z FIRE 

! 


t 
| | 


ot 
i 
+ 





F j 
‘| malt 1 

tte + 4 4 ' = 
yp tee ae i | 

opty etfs | ' 
alg t! : | 5 : 

Vcohiy antl ’ t Saleiaieale| alee 

bes ; 

Amr iOnntins 
| are de aS Re | Teed gene 
SF 4. r trty 
1 . ' ' 








kt 
; : is { ‘ 
’ ah hl | ! ( re | 
pala oil hs fr atom : 
' tygie ‘ ] | “Tir a aaa » | an 
ft efter re ana) | : 
be rer 4 oo - 
' ago ul 1] : “as wt fy ' oar | 


r(em) 


Figure 29 


74 





a, 


PURPOSE 


APPENDIX F 


COMPUTER PROGRAM 1 


1. To correct counting data and compute the saturated 


activity in counts per minute (A 


CD and Ag). 


eee HOmcolnplueeewumMe approximate mean counting rate and, 


LOU onaeveriithneomune Seandard deviation based 


OMPowoOOrMaredrstribpution. 


PEE lomeomouce tae products & 


2 2 
cpr and Agr 


PROGRAM NOTATION 


D 
T 

POS 

cmc 
ASAT, BSAT 
AMEAN 

W 


SIGMA 


Decayare OMe ier Or ini16 


Counter resolving time 

Distanieen rom sceurce 
Leneninmeescomats “side 1, side 2) 
Saturated activity (side 1, side 2) 
Approximate mean counting rate 

Cor reclLaorn Or toni wei echt 
Standard deviation 


Cadmium correction factor 


i) 





APPENDIX F 
COMPUTER PROGRAM 1 


he 


VC aXeVol2 eC CeVY,2Z2 
( 


DISTANCE", 125 5!SAT ACTIVITY ' yT56y "MEANT, 
s 

: 

D 

0 

i 


-— = GN Q ef © ef) enw 
— el Ur py OOO eine SD 
AUS OL et ole He ewe NO) 
R- OF OD J OH ANDO DIMM NO 


REAL*8(A-H,O-2Z) 


mauve NOOO woe <t ll OF Il il 
we w = © eOOOfZ Of e Sere 
CH XM HOMNA YY <P XK MEN It 


am 

om 3t 

NO 

x f 

Cj~ 

ja 

— 

~~ O ui 

sO «xO 

coal tu) 

e Cm 

ir j >- 

MN a 

a) > 

< Oo Cs | 

SD on i(~ 

~ e —~O. 

Ne) uf’ Ax 

nad MN > LL 

<f e Wit 

=. if\ o Co~ 

Oa) ~< ‘a — 

met QO ms! ent 

Ui ~ ~ —, om 

o>< G > —~>x< 

<oO as > ~< + 

<[ © e ~ + © 

Ut .O ry oe) Cyt 

= oi Gi 2) j— 

<{ e QO ~ ~— OO. 

© WU7y ~ MN OX 
Or <li. OM > P- IO) 
«OO WM *©OO & « ~ *O I 
NCO UL KOON? < ee 
~ Ou eQNOMe o@ > OrcIOc) 
—~™ @ @ —- Ol WO Qt QO 
kr Aj O WOQIaiflLn * PCIOO e 
<i“ mt Oe KO Wem mliit, Oo ert 
U1) om St Cina eiGlcy DO OND Aw~ 
COOx LL FF Cee © 6()) 6a 


POT Mmm OOO mm ODO DPD KOON 
Ree ST OR OOD) Pomel em UH ene NC) 
LOBZOWY oF OO™ Fl thi MPreraNnanOouo 
Mm <lOe Owe He RR ODT Ow Ow OY OO 
SSW LL ee he le hee SE AE OE NROCO 
weds Wisi<t celta wee <T HOF It il 
bee Oh RK SONNE Rr OOD OSR KSEE 
PWD ARTO LLL OT CK YN S [I <T 


0 
— 
e 
tN 
GQ! 
a) 
~ — 
~< 0 
ZO oot 
~ e 
Oo iy 
ome OG 
° cS 
Tuy ~ 
=z! x< 
oe xO 
-— & ~ 
MX O 
~O _ 
—e * 
<I:0 Lr 
Ls end NS 
= °@ QO 
BP; <A ~ 
C oI ~< 
O wd (2) am.O 
eOOct ~ ~ a 
NOOMKOO<TN 
NOUN © O00 * e 


Mm + EMO & 

YOR Oh mes ONE 

eT Were Her os 

wer ee TU cD Ld <x ee 
UtekeS qtr akekea 
RL Coa YN OO 


SCOMU POOMOWO™WOWNeTES IOKOUN IDYLOOWOVO— WONKETTE CO WU ~YOCOras! 


sli OrFOOmit atl COOOTE 
re 


ot N OO 


rT @ O mM {= 
uN 


~ MOO 
OO 
uy 


76 





APPENDIX G 


COMPUTER PROGRAM 2 


PURPOSE 
ime LO Compute the area under curved portions of the 
2 2 y 4 . 
Civics Aopr ; Agr ; Aap? and Agr Versus using 


See Orn Smivblleces 

ew LOmecOmMpuue area Under Straight line portions of 
these curves analytically. 

Poe LOPc@mpULCe DOudLmareasvand, from this, Man and 
Try using Equations Gis taearce (ilo). 


4. To compute Tp by correcting T; for additional age 


N 
required to reach thermal energies. 


a To compute dati veven leneth by =solving 


= a /\% 
Lin = (M ee) 


PROGRAM NOTATION 


2 2 y 
y (1) ena values of Aap? » Agr’, Aapr and 
Agr 
H Abscissa interval lengths in cm, used to 
Obtain ycl) valves. 
N Un em oils cits cee mwere ve: IS. i, 
SD) Simpson area approximations 
B(I) Relaxation Length values 
: 2 } y 
ZO) Ordinate values of A,pr’, Agr’, Aap? 


and Aer’ corresponding to abscissa values 
of r at which distributions become exponen- 
alee, 


ch 





A(T) 
AA(T) 
MIG(I) 
AGE(T) 


areal 


Area under exponential portion of curves. 
Omer oemuncder securves. AAC L)= SCJ) + A(1). 
Thermai Migration area values. 

Values of neutron age to Indium Resonance. 


Dimus ron lengths. 


78 


‘ie | Syed ae . a 
cates oe tarot 

ol tsa i® feurgadf? 

tual Th gegld¥ 


gv ; ww ewiyiq 





a 
a Ww 
(14 


: 7 


APPENDIX G 
COMPUTER PROGRAM 2 


Y 
<I 
ee 
oC ~ ~ 
<— ~ ~ i’) —_ 
= mal — <) ome | 
© mr -~+ nN Us + 
Lu CN eee ee 
wy = —- ~n — 
<_{ «<L > 6 a 
ae <t + e+ FT ~ 
er, © cm ° kt eo N -_~ 
mm > = = => ny ~ 
Lig ON Cy. =t2). UC) VY) Be ~ sO 
= + oe Li % = ~ mw 
Ps «<I wd = vn >) BO ® 
<f o (a i it {> rm In 
— a qu ae N LL ° uy oN 
2 o . C F ° (2 LO »e © 
— N Oo = Oo = ON © - 
Co ~~ + Fe © + n 2 cc *~* © 
mM SN a ke oO “NM G&G ~ x HW 
tH - 2 > > Some = ~ =_~ FT 
(1) = > —m— W ~ OO FO m_~ Il =~ os“ a_~ ke 
22 oo NM «KK WZ OW et WN QO <« On Mm © 
mo CGY K ee + © Ue e il .™;“ a 4 
~ OO rd a ee oy O ani tt + << x ~o + a 
> | ~ got > CQ. WY) ~— st Ge oa ~_ — ° ao 
Clo os ™ er ~ ° - de mam my = me 
me - OQ OC ow Oo — ties ~ —- O ~ a —- © 
-r N NI + 2 i" + £2 Ze bent <n e AIKO Utd ome 
ret OY - nm YY ib - a YW JIN <I- —ee x ——- <—[— 2 
<O — ~~ em — ~ tote ~~ e Om =~ —_~ om eS ~™ Cit t~—ith 
Ses Sn 2 CO = x HMY) > O —! Z~HEWMA COW~O mew — oN ert a et S| 
hm LS > iusT = ~+ eK Dm ee FRA! © «> e~K~O~ © DY O om wT (NI fe pede ete OO OR rte 
Ic Oa OSNO 2AAMOUN er omer O2Z+0 2 +UGetU A) met em ewe OU ow OU a Ol 
mM (Oe <t eOi eT O CO) ef Seah Cc eNO OC) Ofer WM LL OEM tert Net NTO NOs eH 
Der OT) Ora ete ONOIO M30 eb tlmesirt-d ONTO OBO eb is<Tei™N er OO We Hed er iu Mt Orr i: Iles efi 
fe D> een I] OO o})/ 4) e HIN Ome Le YP ee CO oi; © WP ed ee ee tee OD ee Lt eg tt Oe Dt IE Oe Dt 2 eS 
LOGWNGMAKRLY Ol OF YWTerZOMNMRrRO OCOUMO YNDHraenreWwn~ers Hho eben | Re Ree we 
VAeON wat N I Wowend Pu WUK<teaeQ wach I ME meh RL ee ee HL TOO Tet Cm i et © I LT 


SIS DOSE *ONONMH 1 OMS THE OCOD HQNANMY | AMR SR OD Or Rr On O era Omer TRO ST A 
COS DOCT LAAN TZOTDIMYZOODALNAST MNT SSOP IDHS AL ZOAIL NH ety ZAR NON YZ OWN SZ rie ZOO 
2 SW wo I! WOOO NOD CYOOZOWO!W 1nDODD NOD ZOO Zwoo~coooatec=eecoceOc ooo} zZeoOr] 
My OTOCXWLZMNOMNEZOIMNSLOMOAYCU TZAWNGNNEOMMN SLO YCUOMFLOOTOT FSLOOSLSULOO KH ZU OMT 


= rt - = 
> © © © OGM a _! ed AOD xe TAN © mr C) eAtC) moO 
Qa =~ ify XY) f- OO. 4 ry O MOA, © Oe Qoo@ Oo 
ae og = ot COUN CVC) rer 
© © © 
W W OQ 


eo 





ReEe NDT Xx 


COMPUTER PROGRAM 3 


PURPOSE 


te 


2 


To determine slope of curve &nA V opto ein che 


tn 
range 18.4 < r < 22 cm. by the method of least 
squares. 


To compute Ton Peon vanes re lat 1. On 


il 


Slope = - We 


PROGRAM NOTATION 


Ee) 


G(I) 


X(T) 


a) 


avi) 


SIMQ 


Relaxation length values 


Values of neutron age and thermal migration 
area 


Distances from the source 
DOmiGlons, CO Equacion 
Ce7r/d 


2 
ie 


De 


Natural logarithm of Z values 
Yo= Ln 24 = ay a Ax X 


suoeoutine fOr SOlVing system of equations 
CAs: 
C A = B 
if if 
; oa 4 “1 ae 
7 aes Lae 
sey “A gE 8 | 2 yey “271 


Slope Of Straight line solution equal 
to the negative reciprocal of Wty 
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APPENDIX I 


RESULTS OF COMPUTER PROGRAM 1 


TABLE 8 


PURE WATER; INDIUM RESONANCE ACTIVATION 





r(cm) Aap (CPM) Mean Count Sev 
emeco 298.0691 O55 elas rele 
9.52 204.1689 oie 3079 12.5454 
On 9 155.0969 2 en 1066753 
re Oah eee oll 89.9505 9.4842 
13.34 omy Or te S 4.3954 
he Cal Cy eos 50.4245 eee eae 
oie 50.9003 BO Soy 622767 
17.14 42,2404 B2.(2L0 Dele oS 
ol 34.2155 25.6226 5.0618 
'PABLE 9 
PURE Wath: Leta AC VATION 
r(cm) Ag (CPM) Mean Count oe 
Seo 4HU59.2573 3602. 2972 60.5252 
Oe 52 3549.0267 Boe. cine 54.0207 
10.79 Afi se il 2259.8505 Wes 246 
ero 1930.5774 leeds ce22 Ono 2 
13.34 ey eho il 1181.5408 Seles 30 
eo LOCO, 6 202 Ole alec) SO Vea. 
15.88 852.3489 692. 3284 2G t 
ieee 642.0987 eo Ll sy Pee 
eed. 5 eS 3 449.0388 Pale US 
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iO 
3 


TABLE 10 


CONCENTRATION 1.16 g/1; INDIUM RESONANCE ACTIVATION 


3 

i Standard 

r(cm) Aap (CPM) Mean Count ee 
8,26 503.1124 RTOs 19.4766 
9.52 pecre eae is 240.9674 ese ceoL 
105.7 9 2 oe Oe IeEC OOM: S 13,2665 
ez Oy, Oe oO ly le ae ese eek 
Ae ee oO 2 Oe Se ie seue > 
eo leer adel 90.3156 Oa 5) 
15.88 74.4533 57.2474 7.5662 
ieee GNNG303 em ions 8.5196 
ore ees 36.5259 6.0436 

AAG eels 

HBO. CONCENTRATION 1.16 g/1; TOTAL ACTIVATION 
Standard 

r(cm) A.(CPM) Mean Count ee ea ae 
Cara 6261. 3213 le Oe eu Oe 
Oma Hee relel: 3446 .4026 58.7061 
LOR 79 3113.6634 2596. 8206 Oe OS 30 
ZO, POO too 1854 .0849 43.0591 
13.34 Poole 5 o2° 1400.5575 Bae 
ee oi iis. oO 1080.4300 32,8698 
15.88 986.9352 CH ee fae 29.0254 
eet 768.5511 639.9204 25.2966 
Moe 635.0632 Ber. poo eee lay 
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H_BO 
3 


TAbee. 12 


3 CONCENTRATION 8.96 g/1; INDIUM RESONANCE ACTIVATION 
r(cm) A op (CPM) Mean Count oo... 
8.26 409.9467 314.5537 if 35 0 
9.52 oie ia PMs 2628 14.6035 
Oey 229.5099 on embed 12.9233 
eas ogee eee 123.1536 11.0975 
13.34 12 BOs MOB ob iV Onere t 
Poa. 105.7302 81. 3866 S20 ealS 
Sos Geo 51,4156 fe a0 
ieee Ser one 62.9505 7.9345 
eer al 44.4031 33.5245 5.7900 


AAG es 


HBO. CONCENTRATION 8.96 g/1; TOTAL ACTIVATION 


r(em) A,(CPM) Mean Count = Sieur 
8.26 Pleo boa. 2326. 7305 W236 1 
oe LRP OHS) 11 steely 37.5961 
10.7.9 eee 5005 986.2288 sie 4042 
es On7 Salas © enks 725) 0108 26.9333 
13.34 639.1190 530.0796 23.0235 
een 487.743 405.6754 2 One 
See © 373-6443 SiO. aie 74 infeo2e4 
17.14 2 Oe Slt 226.8509 15.0616 
ee 269.3002 220 70e22 14.8334 


By 





TABLE 14 


aS CONCENTRATION 27.81 e/1; INDIUM RESONANCE ACTIVATION 

Standard 

r(em) Map (CPM) Mean Count ean 
8.26 285.2747 218.7960 ieee 8 
OR De 236, 3901 183.6949 13.5534 
AOE FES Oreo lee iL espa he Wie 1 
207 a Sree Ba NO 3234736 Oz 2 
13.34 146.9174 Oa) sie 3 Ores ec 
ar eGo 89.3376 66.5402 Be oye |. 
15.88 69.3343 palo es ole 
17,14 93.9194 (omer eo oe 
18.41 42.0807 Be aloe 5.7049 

ect) eel 

HBO. CONCENTRATION 27.81 g/1; TOTAL ACTIVATION 
Standard 

r(em) Ag (CPM) Mean Count eee on 
8.26 a3 3 3 0 PU bocs 33.4450 
Ore 855.6882 eo ae oA eo ln 
Vor. 78 658.1651 Bison lay 2 1 PB 0C50 
207 BOG a 4HO7.9465 20.6868 
ts. 34 360.5809 306.2933 gaa Oe 
eon 399.0305 Be ecole], Wee cece 
15.88 206.7903 76m oes 9 sec oo | 
lieve 186.5858 158.2384 hee 1 93 
gers cat 1297, 9547 Lies TS / Oe 50s 
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TABLE 16 


HBO. CONCENTRATION 46.65 g/1; INDIUM RESONANCE ACTIVATION 
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a 3 
Standard 
r(em) Aap (CPM) Mean Count es 
8.26 eon ies Oe TARO. 17.2543 
O52 282.4397 Ole e No 1S 76 
iO <7 9 AO elo Zo 153.5804 een, 
2. OH oe. OO ALO 25 10 lO 26524 
13.34 LOSS 517 113.8245 10.6688 
aol) O0rose2 | 67.9496 8.2432 
15.88 95.1983 (ibe Ss 8.4768 
1 fad 54.4682 40,6000 Onsale 
Hee oS hy eo fee 7.5946 
TPABEE. 17 
HBO. CONCENTRATION 16.65 g/1; TOTAL ACTIVATION 
standard 
r(cm) Ag (CPM) Mean Count ee aa 
Oneo OGo. 15> / 835.3732 2G OG 
OF a Onion 3k 2 Se 23.2922 
Ono ers © 5) Bn OOlLg LOS4 367 
Oy. 364.3542 202 720 ey eee 1 
13.34 310.0154 251.2099 15.8496 
14.61 22 Cas 5 Gone Lee 13.6458 
15.88 199.8872 159.1900 1236472 
Nae ieee o 2 onl Ole Wiel Gy Selo 
16.4) Bo. 50M 73.0562 Cea 





WUE O 


HBO CONCENTRATION 70.09 g/1; INDIUM RESONANCE ACTIVATION 


is 
Standard 
r(em) Aap (CPM) Mean Count ee 
cece 379.0630 285.6253 16.9004 
9.52 POV TONS ne ome ee oe 
Ono 230.4439 165.9322 orem 
mann 160.5295 120.9588 10.9981 
13.34 106. 3356 79.5638 8.9198 
eno) 95.0191 69.8891 8.3599 
eo ee 55.6762 Paes? 
ane Slesicie 38.2288 pie 2° 
iS a 54,6934 a One Se eal 
TABLE 19 
HBO. CONCENTRATION 70.09 Ce ieatoli ACTIVATION 
Standard 
r(em) Ag (CPM) Mean Counts eae 
8226 824.4368 683.2489 DE Be 
Gr 52 558.2886 462.7575 ae > aleay 
10.79 496.2032 aoe eo 16.0206 
a oeaOnr 316.5702 261.6614 16.1759 
13.34 PE Gnomes 189.6963 Seen 
een! ee olen SO plo / leat 2 
15.88 137.7566 ae 10.6851 
a ade 93.6648 79.3185 CO Ou: 
Poe FE OGL 60.5687 Veeco 
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TABLE 20 


HBO, CONCENTRATION 80.09 g/1; INDIUM RESONANCE ACTIVATION 


r(cm) Acp(CPM) Mean Count Sea 
8.26 414.8513 307.2195 oer 
9.52 256.6541 196.2806 14.0100 
Ona Pere iel 177.5645 13.3253 
ee Og 161.0497 122.6576 Oneal 
ieee 113.0485 87.4346 9.3506 
oul 99.8647 71.9423 8.8285 
15.88 90.3185 68.9324 8. 3026 
17.14 60.4457 45.0622 6.7128 
eo eat 35.4877 27.8769 5.2799 
TABLE 21 
HBO. CONCENTRATION 80.09 g/1; TOTAL ACTIVATION 
r(em) Ag (CPM) Mean Count See pareja 
8.26 743.8847 614.4305 24.7877 
9.52 No, 7427 419.9927 Ag os7 
10.79 SOC nse SOS) 5 oe S5) ogee 
Tao 302.0725 261.4078 16.1681 
13.34 Pale oe 175.3501 ee aC 
14.61 149.1620 128.0910 11.3177 
15.88 ie eae 105. 3823 10.2656 
17.14 116.5993 99.1096 9.9554 
18.41 70.7709 53.4546 ees dal 3 
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